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Abstract 

We regularize in a continuous manner the path integral of QED by construction of 
a non-local version of its action by means of a regularized form of Dirac's 5 functions. 
Since the action and the measure are both invariant under the gauge group, this regu- 
larization scheme is intrinsically non-perturbative. Despite the fact that the non-local 
action converges formally to the local one as the cutoff goes to infinity, the regularized 
theory keeps trace of the non-locality through the appearance of a quadratic divergence 
in the transverse part of the polarization operator. This term which is uniquely defined 
by the choice of the cutoff functions can be removed by a redefinition of the regularized 
action. We notice that as for chiral fermions on the lattice, there is an obstruction to 
construct a continuous and non ambiguous regularization in four dimensions. With the 
help of the regularized equations of motion, we calculate the one particle irreducible 
functions which are known to be divergent by naive power counting at the one loop 
order. 
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1 Introduction 



As it is well known gauge invariant quantum field theory is a powerful tool to describe 
fundamental interactions in physics. However these theories suffer from ultraviolet (UV) 
divergences if we impose the locality of the interactions. Over the past decade since the 
birth of QED great efforts have been made in order to regularize these theories in a gauge 
invariant manner in order to preserve the Ward-Takahashi (WT) identities, since it appears 
hopeless to cure this problem by the construction of non local interactions [|TJ. This is a 
necessity dictated by the renormalization program. There are essentially two categories of 
gauge invariant regularization schemes, a discrete one and a continuous one. In the first 
category the path integral of the theory is discretized on a lattice. This allows to make non- 
perturbative, but mostly non analytical predictions. On the other hand the second category 
of continuous gauge invariant regularization schemes is intrinsically perturbative, in the 
sense that in momentum space these schemes basically regularize the Feynman integrals or 
the lowest order Green's functions in configuration space. This is the case for Pauli-Villars 
regularization ||, dimensional regularization || |J, Schwinger's proper time methods inspired 
regularization |J, ( function regularization ||, differential regularization |7j and Fujikawa's 
type of regularization |^. |9| . 

If one wants to study the non-perturbative properties or the different phases of a gauge 
theory as a function of the scale, a general method is the Wilson-Kadanoff Jl0| approach. In 
this case it is more sensible to work, in four dimensions, with a continuous regularized and 
gauge invariant version of the path integral of the theory. In recent years, some progress 
have been made in this direction. For instance, gauge invariant regularization of chiral gauge 
theories can be achieved by a generalization of the Pauli-Villars method [11 , || i.e. by the 
addition of a finite or infinite set of non physical regulator fields at the lagrangian level. 
In the same spirit a generalization of Schwinger's proper time methods, the operator cutoff 



regularization ||12|| , allows to calculate effective gauge invariant regularized action at the one 
loop level. 

In the present work we follow Wilson' s idea, but we regularize non-perturbatively the 
path integral of QED in a continuous manner. This is done by a suited smearing of the point- 
like interactions of the fields by cutoff functions which are a regularized form of Dirac's 5 
function. We choose the cutoff functions to act only in the (UV) domain and not in the 
infrared (IR) one. In this sense, for finite (UV) cutoff our action is non-local, but converges 
asymptotically to the local one as the cutoff goes to infinity. We should notice that an earlier 
attempt was made in constructing a non-local regularization of gauge invariant theory ||13| . 
In this regularization scheme only the action is invariant under a non-local transformation^ 

lr rhis transformation is not the representation of a group. 
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which plays the role of the gauge transformation. The invariance of the path integral is 
maintained by adding a term in the action which cancels the contribution of the transformed 
measure. Since this term must be worked out order by order in perturbation theory, this 
method belongs to the category of perturbative regularization schemes. 

In our case, both the cutoff action and the measure are invariant under the gauge group, 
therefore there is no need to define the regularized partition function perturbatively through 
the regularization of the Feynman diagrams. As a result the path integral is regularized 
as a whole in a gauge invariant manner. This feature gives us the opportunity to apply 
directly non-perturbative techniques to the regularized partition function in order to derive 
the physical properties of the theory. For instance, in addition to numerical simulations, we 
can study the variation of the path integral under an infinitesimal variation of the cutoff 
scale in order to derive Polchinski's [0] like exact renormalization group equation. We can 
also derive the regularized equations of motion, i.e. the regularized Dyson-Schwinger (DS) 
equations, directly from the regularized partition function of the theory. Since this scheme 
keeps track of the dimension of any regularized integral through the logarithm or the powers 
of the cutoff, the domain of validity of the regularized (DS) equations is not only restricted 
to the (UV) region, but can be extended to all range of the cutoff scale in contradistinction 
to dimensional regularization ||. Thus this regularization scheme can give an insight into 
non-perturbatives properties of the theory, even in the context of perturbative calculations 
where the (DS) equations appear as a resummation of perturbative series. For example this 
scheme is suitable for the calculation of the vacuum expectation value of operators like the 
condensates. Moreover, as this method is continuous and works strictly in four dimensions, 
it can be applied to the treatment of the axial anomaly which is known to connect both the 
(UV) and the (IR) regions. 

In order to show that the partition function of (QED) is indeed finite, we derive first 
the regularized Dyson-Schwinger (DS) equations which are the basic ingredients for the 
skeleton expansion. The class of regularization is imposed by the choice of a cutoff function. 
We choose to mimic Schwinger's proper time regularization. Then, despite the fact that 
our regularization is non-perturbative, we deduce the regularized form of the relevant one 
particle irreducible (1PI) functions which are known to be potentially divergent by power 
counting^ perturbatively from the regularized (DS) equations. We show that we recover the 
known result for the mass operator. This is also true for the vertex function if we assume 
that the electron lines are on mass-shell. Since the regularization of the path integral is gauge 
invariant by construction, the polarization operator is transverse and regularized. The finite 
part of the polarization operator has the standard form, the only discrepancy occurs in its 

2 The direct study of the regularized partition function with group renormalization techniques will be 
done in a forthcoming paper. 
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divergent part where a quadratic divergence appears. This divergent term can be removed 
by adding to the lagrangian a gauge invariant counterterm S\ which is quadratic in the 
photon field. The paper is organized as follows. In section 2 we define the regularized gauge 
invariant action in configuration and in momentum space. The regularized form of the (DS) 
equations and of the (WT) identity are given too. We derive the relevant (1PI) functions 
from the equations of motion in the next sections. Section 3 is devoted to the polarization 
operator and to the 2n-photon amplitudes. In addition to this section the role of Sa is 
clarified. We calculate the electron mass operator and the vertex function in section 4. 

2 The regularized gauge invariant action, the (DS) 
equations and the (WT) identity 

In order to regularize the action of QED in four dimensions we construct a non-local cutoff 
action which converges asymptotically to the standard non-regularized one as the cutoff 
tends to infinity. This can be done by a suited smearing of the point-like interactions of the 
fields by the scalar cutoff functions 

Pi (x,y) = J dke~ ik ^ Pi {k,A), (2.1) 

which are symmetrical and are regularized forms of Dirac's 5 function, i.e. 

lim Pi (k,A) = 1, (2.2) 
lim pi(x,y) = 5{x-y). (2.3) 

Here and in the following we work in Minkowski space. We choose the signature of the metric 
to be (1, —1, —1, —1) and the notation dx = d A x and dk = d A k/ {2ir) A . 

In all analytically methods used to calculate or to define the Greens's functions of a 
given theory the inverse of the energy kinetic terms appear through the propagators. This 
is the case for non-perturbative methods based on exact renormalization group equation or 
Dyson-Schwinger expansion based on the use of the equations of motion of the theory and 
or pure perturbative methods based on Feynman diagrams expansion. It follows that the 
regularized action S# eg of QED must obey three conditions. 

(1) SR e g must be gauge invariant. 

(2) In any expression composed of product of (1PI) functions the inverse of the cutoff 
functions ( j2.ip associated to the kinetic terms which occur in the propagators must not 
cancel the cutoff functions associated to the vertices. 

(3) When the cutoff A goes to infinity Sn eg must converge to the local non-regularized 
action of QED. 
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Since in QED the photon has no self-interaction we need only to regularize the fermionic 
part of the action. A solution to these constraints is to regularize separately in a gauge 
invariant manner the kinetic, the interaction and the mass terms. The minimal solution is 
to keep the mass term unregularized. Then the building blocks needed to construct Sn eg are 
the smeared gauge field 

A"(x) = J dy p 3 (x,y)A^(y), (2.4) 
the functional of the smeared gauge field 

L{z,z') = J dyC,(z,z',y)A^y) (2.5) 

and the functionals 

*i(x) = |dy Pi (x,y)e ieL ^V(y) (2.6) 

of the smeared fermion and gauge fields which remembers Schwinger's point splitting ||. 
Likewise in lattice gauge theory the functional L(z, z') plays the role of a link. The form of 
the C-number function is imposed by the necessity of gauge covariance. In order that 
the functional ^i(x) transforms covariantly under the gauge transformation 

A^x) - - ifyA(z) 

i)(x) -> e l I dy p3(x ' y)Aiy) ij(x), (2.7) 

the variation of L(z, z') must be 

SL(z,z') = ijdy (p 3 (z,y)-p 3 (z',y))A(y). (2.8) 

This imposes that the divergence of the C-number function verifies 

d^C^zJ.x) = 5(x-z)-8(x-z'). (2.9) 



In addition the choice of is such that 



L(z,z) = 0. (2.10) 



This is the necessary condition to recover the local fundamental fields from ( 2.4j) and ( |2.6j ) 
when the cutoff tends to infinity. In fact due the properties (|2.3|) of the cutoff functions we 
have now, 

lim A M (x) = A"(x) 

A^oo 

lim * s (x) = (2.11) 

A^oo 

For we choose the less (UV) divergent solution 

C^z, z', x) = % [ 'dr% ( e —(— ) - e- ir ^- z,) ) . (2.12) 
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Then the regularized action Sn eg of QED is the sum of the fermionic term SWmion and of 
the pure gauge term Scauge which are defined by 



S 



Fermion 



dxdzdz 4>(z) e ieL{z ' x) 



ep 2 (z,x)[J)L(x.) + (#L(x,z')) )p 2 (x,z / 



pi(z,x)(i(p + e(^L(x,z')))p 1 (x,z') 



3 ieL(x,z') 



— m 



dx ipip 



S, 



Gauge 



dx 



+ ^A,f] + S. 



'A- 



(2.13) 
(2.14) 



The contribution S\ of ( 2.14|) is gauge invariant and will play the role of a counterterm for 
the polarization operator. This term will be discussed in more details in section 3. Since 
L(z,z') is defined by fl2.5|) and ( [2.12j ), the partial derivative $L(x,z') which enters in Q2.13|) 
is a function of x only. This implies that the kinetic and interaction terms of the regularized 
action (|2.13|) have the general form / dx ^jKi^i, where Ki is a functional of the smeared 
gauge field ( |2.4j ). As a result one can readily check that each of the three terms of ( |2.13| ), 
which correspond respectively to the kinetic, the interaction and the mass term, are by 
construction separately gauge invariant. When the cutoff A tends to infinity, it follows from 
the relations ( |2.11| ) that the terms proportional to (j)L(x, z') cancel and that the regularized 
action Sn eg — S\ converges formally to the standard local action of QED. In deriving the 
relations ( [2.1 1|) and the asymptotic form of the regularized action we have assumed that we 
can formally interchange the limit with the integral symbol. This assumption is not exactly 
true0, and will be to the origin of the appearance of a quadratic divergence in the polarization 
operator. If we use the commutation relation 



a ieL(x,z') 



— e e ieL ^ x ' z '^ 



{$L{x,z% 



(2.15) 



we can rewrite ( |2.13|) in a more concise form as 



S 



Fermion 



dxdzdz' ip(z) e ieL(z ' z,) 



Px{Z, Xffl p!{X, Z 



- ep 2 (z, x) L#(x 



m I dx ipip. 



*L(x,z')) Wx,z') 



V(z') 



(2.16) 



From now on, the regularized gauge invariant path integral which we consider in configuration 
space is 

Z{rj fj J) = [ VipVtpVA^ e^ 5 ^' - ^'^ + J dx (fi(x)ii)(x)+i>{x)n(x)+j^(x)A»(x))) ^ 17) 



! Evcn in one dimension this is not in general the case for improper integral. 
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where S(ijj, ip, A) is the sum of ( [2.14 ) and (|2.16|) . The T], fj and J a are respectively the 
sources for the fermion and gauge fields. As before one can notice that the new form ( p. 16 ) 
of the fermionic part of the action is still composed of three terms which are separately gauge 
invariant. The fact that the whole action S(ip, iJj, A) is a regularized action of QED in four 
dimensions will become transparent in momentum space where we give the representation of 



the partition function (|2.17|) and of the Legendre transform of the (1PI) generating functional 

G{r], fj, J) = log Z(r},f),J). (2.18) 
For the partition function we have 

Z(n,fj,J) = J VipVipV A a e i ( s( ^'^' j4) + $ dl (v(kW(k)+$(k)ri(k)+J ll (-k)A> l (k))) ^ (2.19) 

In this expression the action is 

S{ip,^,A) = j dkdpdp' 4>(p)K(k lPl p')il>(p') + S Gauge (2.20) 

Saauge = ~\j dk (&V - k u k u (l - ±)) A»(k)A V (-h) + S A , (2.21) 

and the kernel K is given by 

K(k,p,p') = F(p-fc,A;-p')p?(^)^-^((27r) 4 ) 2 5(p- k)5(p' - k) 

-ejdq p 3 (q)A' l (q)T^q)F( P -k,k-q- p')p 2 (k)p 2 (k - q), (2.22) 

where the transverse matrix T^q) f\ is defined by 

r M (g) = 7m -i% 



(2.23) 



and F(p,q) is the Fourier transform of e teL ( x,z '\ With the help of the recursion formulas of 
appendix B which define F(p, q), we can rewrite the action Q2.20|) in the final form 



dk ip(k) (pi(k)$ — mj ip(k) — e / dpdp' ip(p)A^(p — p')T IJr (p,p , )ip(p') 



~T In 



n=0 



(ie)("+ 2 ) 
(n+2)! 



dkdpdp' ip(p) 



F n+2 {p-k,k-p')p{{k)^ + t(n + 2 



J dq p 3 {q)A"{q)T^q)F n+1 {p -k,k-q- p')p 2 (k)p 2 (k - q) 



Here the matrix T^p, p') is defined in terms of ( |2.23[ ) by 

r„(p, j/) = P3 (p - p') [p 2 (p)p2(p')t,(p - p') + ( P l(p)t> - plW) 



(2.24) 



(2.25) 



'See the properties of this matrix in appendix A. 
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(a) (b) 
Figure 1: The vertex associated to the matrix — eT^p, p') (a) and to its transversal part (b). 

and has the property 

lim r„(p,p') = 7 „. (2-26) 

A — >oo 



Owing to the property ([2.26 ) the first and the second term of fl2.24|) converge to the standard 



fermionic part of the action of QED when the cutoff goes to infinity. As for the third term, 
it vanishes identically due to the properties ( P-6| ) of the F n . This term which describes an 
infinite set of vertices is needed to ensure the gauge invariance of the regularized ampli- 
tudes. Therefore the whole action ( |2.24|) converges formally to the non-regularized action 



of QED. Since the third term of (|2.24|) vanishes asymptotically, we expect that the matrix 
— eY^p.p') ( p.25|) or its transversal part —ep? J {p~v')P2{p)P2{p')^^{p — p') associated to the 
vertices of Fig.|TJ will play an essential role in the representation of amplitude by Feynman 



diagrams. Notice that the matrix r^(p,p') (|2.25|) is unambiguously defined by the values of 



the momentas of the two electrons lines because ^^{p^p') is symmetrical in the variables p 
and p'. 

From the action ( |2.24|) it is readily seen that the free electron and photon propagators 



behave respectively for high k 2 like 1/m and 1/fc 2 . Since at each vertex is associated a cutoff 
function, it follows that the path integral ( |2.19| ) defined by the action ( j2.24 ) is finite. Apart 



from the fact that in euclidian space the (UV) cutoff functions must be rapid decreasing 
functions of the squared momenta and must be defined for all values of the momentum, 
the choice of their form is quite arbitrary. For instance if all the cutoff functions ( |2.1|) are 
identical the action (|2.24|) which is expressed in terms of the fundamental fields does not 
simplify. However when expressed in configuration space in terms of the smeared fields ( |2.4| ) 
and ( |2.6| ) the fermionic part of the regularized action take the compact form 

= /tJi^-m/dxfe (2.27) 
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where iJfi is the smeared covariant derivative defined by 



ip = e#. (2.28) 

From a practical point of view it is better to leave the cutoff functions unconstrained, then 
in the following we consider the general action ( |2.24j ). In order to check that the theory is 
indeed finite, we shall derive in the next sections the regularized (1PI) functions from the 
regularized equations of motion. 

For this purpose we also work out the Legendre transform of ( |2.18| ) in momentum space. 
It is given in terms of the vacuum expectation value of the fields 

m = *(27r) 4 4ty m = -K^Tsm M-k) = -Wjlk) (2-29) 

by 

T(M,A) = -iG(r),fj,J) - fdk(n(k)i/)(k)+^(k)rj(k) + J^(-k)A^(k)y (2.30) 
We recall the conjugated relations for the sources 

V(k) = ~ (2vr) 4 ^ v(k) = (2vr) 4 ^ J,{-k) = - (2tt)^. (2.31) 



From the invariance of the partition function ( 2.19 ) under infinitesimal translation of the 
fermion and gauge fields and under infinitesimal gauge transformation we deduce respectively 
the (DS) equations and the (WT) identities which are now mathematically well-defined 
objects. The regularized form of the (DS) equations is 

M-r) , SS A 
(2tt) 4 5A V (r) 



+ '-M + Sh) (2-32) 



= (/ dkdp 1 K(k,p,p')ip(p') + rj{p)), (2.33) 

where the kernel K(k,p,p') is explicitly and implicitly given by the expression of the action 
( ^.22| ) or ( |2.24j ) respectively. Due to the relations ( ^.31 ) the regularized expression of the 



(WT) in terms of the generating functional of (1PI) functions can be written in momentum 
space as 



(2t) = ep 3 (r) J dk [m^h) + sm^k ~ r)\ - \r^A,{-r). (2.34) 

Except for the additional cutoff function p 3 the above expression is identical to the standard 
one. Thus we expect that all general relations among Green's functions which can be deduced 



from ( 2.34|) are the same as in the case of known regularization schemes. 



3 The polarization operator and the 2n-photon ampli- 
tudes 



In many gauge invariant regularization schemes the calculation of amplitudes relative to 
closed fermion loops is a test of the method. In fact any contribution to these amplitudes 
must be unambiguously finite. We will see that this is indeed the case for the polarization 
operator and that due to explicit gauge invariance of the method this conclusion remains 
valid for the 2n-photon amplitudes. 

In order to calculate the polarization operator we take the functional derivative of the 
equation of motion (|2.32|) with respect to the C- number function A^{— r'). Then using ( p.31|) 



to express J u (—r), the inverse of the full photon propagator flC.3j ) is given in terms of the 
polarization operator ir^ by 

( 2n ) 4 6A»(-%A»(r) = 7r ^( r ? r ') - 5 ( r -r')(r 2 g^-r^r v (l-\)) 

l S i SS A \ 
SAv(-r') \5A v {r)' 



+ (2vr) 4 ^y(^) + 0(A) (3.1) 



r,r') = (2vr) 4 - f -^ / dkdpdp' (^(p)^g^(p'))- (3-2) 



Here O(A) represents terms which are vanishing when the sources are switched off. We 
will now calculate the right hand-side of ( |3.2| ) perturbatively at the one loop level. The 
comparison of this result with the known expression of the regularized form of will 
then constrain the form of the operator S\. Since in the action (|2.24j ) each monomial is 



proportional to the product of n gauge fielcf] and is of order e n , at order e 2 , the operator 
ip(5K I '8A U )%1) is at most linear in the gauge field. Hence using (|2.20|) and fl2.24|) its vacuum 



expectation value can be expressed in terms of the full electron propagator 

S(p',p) = - i (2nY w ^ M , (3.3) 

which is defined in presence of external sources, and only in terms of the vacuum expectation 
value of the gauge field. 

Knowing that the functional derivative of S with respect to A^{— r') is at least of order 
e, we obtain for 

7v(r,r') = ^Tr^eJdpT^p-rY-^lA 

- \e 2 J dkdpdp' Ukpsir') (ja%)(p -k + r',k- p') 



-^{p-Kk-pi + ri))p\m + P2 (k)p 2 (k + r , )p 3 (r')T l ,(r') 



<5Fi 



vhip- k,k + r' -p') + (r' -r, p, <-> v) 



SA v (r) 



S{p',p)\ (3.4) 



5 See the expression of F in appendix B. 
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If r< 2 > is the inverse of the full electron propagator, we can express in terms of the three 
point function ( |C.5|) , and using ( p.4j) and (p.5|) we get at order e 



(3.5) 



Here r M (p',p) is defined in ( ^.25[ ) and S(p) is the free electron propagator whose expression 
is readily deduced from the first term of ( |2 .24 ) and is given by 



S(p) 



(3.6) 



Due to the properties (|2.26|) and ( P-6| ) we can notice that the expression ( |3.4|) converges 
formally to the non regularized form of the polarization operator when the cutoff A goes to 
infinity. 

Now if we substitute in (|3.4j) 5S/5A 11 by its expression ( |3.5|) , we can see that the polar- 
ization operator 



(3.7) 



has the following structure 



V,(r) = Tr[r^r u A + (r^T u (r)B 1 + r^(r)5 2 ) + T^r) ® Y v {r)C] . (3.8) 



Since the theory is gauge invariant by construction, the (WT) ( |2.34j) tells us that the polar- 
ization operator must be transverse. Owing to the transversality property ( |A.2|) of the T M 
matrices, this implies that the matrices A, Bi and B 2 of ( [3.8D must be zero if we multiply 
both sides of ( |3.8j ) successively by r^r v , and then by r M and r u . This is indeed the case. A 
straightforward calculation shows that B\ and B 2 are both proportional to the integral 

(3.9) 



(3.10) 



dp P2{p) {P2{p + r)S(p) - p 2 {p - r)S(p - r)) , 
and that the matrix A is proportional to the integral 

dp (pi(p)^ — m ) S(p — r) — (p\(p + T)(J» + f) — rnj S(p) 



Since for high squared momenta the asymptotic form of the free electron propagator ( |3.6| ) 
is 1/m the two integrals (|3.9|) and (|3.10|) are regularized and vanish identically after a shift 
of variable. However there is a great difference between the matrix A and the matrices B\ 
and i?2- Whereas in the massless case one can always choose the cutoff functions in order 
to regularize the matrices Bi and B 2 given by ( |3.9|) , this is not true for A where now terms 
like pf(p)pi 2 (p — r) appear in the integral (|3.10j ). In contradistinction to the massive case, 
it seems very difficult to construct a continuous non-perturbative regularization scheme for 
massless fermions which is also mathematically clean. This is reminiscent of the Nielsen- 
Ninomiya Theorem \TE\ which implies that a chiral invariant regularization is still lacking in 
the fermionic sector of lattice gauge theory. 
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k + r 

Figure 2: The polarization operator. 



We thus obtain for the regularized form of the polarization operator (372) 



( r ) = l j£fpl( r ) Tr I dk 



p\(k)p\(k + r)T^r)S(k + r)T u {r)S{k) 



(3.11) 



The polarization operator can be represented by the diagram Fig.[| where all the vertices are 
transverse Fig{l]b. In the present work we calculate the relevant (1PI) functions with cutoff 
functions of the form 

k 2 

p{k) = eA% (3.12) 



A being the (UV) cutoff. In order to simplify the calculation we choose, 

p\{k) = p 2 {k), p 2 {k) = p 3 {k) = p{k). 



(3.13) 



In this case we can easily show that 

fdkf(k,A)S(k) = fdkf(k,A)p- 1 (k) ^! +m + O(f) (3.14) 

J J k 2 (l+-^)-m 2 

provided that \f(k, A)| behaves like p 2 (k)/k for high k 2 . This implies that we can always 
make the replacement 

S( k )^(l-^)p^(k)-^ +m 



k 2 —m 2 



(3.15) 



in the integrals which we encounter in the calculation of the relevant (1PI) functions without 
changing the result of the integration. 



If we take the trace after the substitution ( |3.15|) and use the Feynman parametrization 
for the product of propagators the expression (|3.11|) can be written as 



Aip\r)^{l-2^)f Q dxjdk 



[—k 2 —2krx+m 2 —r 2 x] 2 



' 3k^_ 3 (k+r) 2 

e 2 a 2 e 2 a 2 



(At) r^r u - r kr{r^k v + r y k^) + r k^k v - \(r g^ v - r^r u ){k + kr) 

k 2 {k+r) 2 ^ 



+eA 2 e a 2 m r (r g^ v — r ^r v ) 



(3.16) 
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The asymptotic form of the integrals which contribute to the right hand-side of ( |3.16| ) are 
given in appendix D. With the results (p.7|) and (|D.8| ) and if we introduce an arbitrary unit 
of mass /i, ( p. 16 ) can be rewritten 



2 / rfx x(l - x) log "'-^d-*) k (3.17) 



o 



a 



being the fine structure constant. We see from the expression Q3.17 ) that the polarization 



operator is transverse as expected and has the standard form |T7|, except for the quadratic 
divergent term which is proportional to 1/r 2 . As noticed in Section 2 the origin of this term 
is due to the fact that in configuration space the regularized action SR eg — S\ (|2.13|) and 
( |2.14 ) converges exactly to the non-regularized action of QED if one makes the assumption 



that taking the limit commutes with integration. Stated differently, this discrepancy is due 
to the non-uniform convergence of the integral in momentum space. The new divergent term 
of ( |3.17 ) can be removed by adding the counterterm 



Sa = £ Jdkpt(k)^{k% v - k^K) [f -m 2 (|-21og|)]^(A:)A l '(-A ; ) (3.18) 

to the expression of the regularized action S^g in momentum space. The counterterm Sa 
can be viewed as associated to a new bare interaction 

S c = c^J dk pl(k)±{k 2 g^ - k,k u )A»(k)A»(-k), (3.19) 

where the coupling constant c has the dimension of the square of a mass. This interac- 
tion which is non-local in configuration space can be absorbed in the inverse of the photon 
propagator ( |3.1| ). Since at order a n the photon propagator ( |C.3| ) will now contain a term 
proportional to c n / {k 2 )^ n+1 \ the contributions of this interaction to the potentially divergent 
(1PI) functions are all (UV) finite. This implies that the new coupling constant c, which 
can give masses to the photon, is not renormalized by higher order corrections. It follows 
that the counterterm Sa is uniquely defined by the choice of the renormalized form of the 
coupling constant c. In order to keep the photon massless as in standard QED and to avoid 
any non-local interaction in the action, we restrict ourselves to the form ( |3.18| ) of the coun- 
terterm, which means that the renormalized form of the coupling constant c is tuned to zero. 
This gives the known result for the polarization operator (|3.2j ). Another choice for the bare 
coupling constant c, that is to say for the counterterm Sa, seems to open the possibility of 
dynamical mass generation for the photon through the non-local interaction (|3.19|) . 

Now we show that the expressions for the 2n-photon amplitudes are standard. We begin 
by the calculation of the polarization tensor. If we use the equation of motion ( |2.32| ) and 
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the definitions ( |2.31| ), since Sa is quadratic in the gauge field, the (1PI) function relative to 
the four photons amplitudes is given by 

8Af i {r 4 )SA a (rt) 5 A" (r 2 )5A^(r 1 ) = $Af> l (r 4 )SA a {r 3 )5A" (r 2 ) Tr J dkdpdp' SA „^ (k, p, p , jj) 

S(p',p)e G + 0(A), (3.20) 

where the differential operator 8K/8A* 1 with respect to the gauge field sources J is easily 
obtained from ( ggg ), ( ggg ), ([K3Q and (Q. 

In the last expression it is understood that 5KJ 8 A acts first on the product Se G and then 
that we act successively with the functionals derivatives 5/5A. Knowing that the action of 
8/8 J or 8/8 A on S ( |3.3| ) gives rise to terms which are at least of order and that at the 
one loop order terms of the order C?(] v ) are neglected the polarization tensor 

riSa^(n,r 2 ,r3,r4)«J(ri+r 2 +r3 + r 4 ) = ((2?r) 4 ) 3 ^ (r4)(5Aa(r ^ (r2) ^ (T , i) (3.21) 



can be written as 

r( Sap( r t,r 2 ,r 3 ,r 4 ) = Tr 



TfATrtro&rpiA + (r v2 r a3 r / 34T^(r 1 )Bi + three similar) 
+ (r^r^T^r^ g) T u (r 2 )Ci + five similar) 
+ (?>JV(ri) ® Y v {r 2 ) <S> T a (r 3 )Di + three similar) 



+ r M (n) <g> r„(r 2 ) ® r a (r 3 ) ® r>(r 4 )£ 



(3.22) 



Here A, Sj, Cj, -Dj and E are scalar functions of the four momenta ri, r 2 , r 3 and r 4 and 
tensorial expressions with respect to the spin indices. Since the regularized polarization 
tensor is gauge invariant by construction, it follows that A is identical to zero if both sides 
of ( P-22 ) are multiplied by the product r^r^r^r^. Similarly we can show that only the 
scalar function E in (|3.22|) is non vanishing. The structure of T^ a/3 being now established, 
we keep only the term proportional to in 8K/8A. Since this term is not a differential 
operator the calculation of the right member of ( |3.21| ) is indeed reduced to the evaluation of 
8/8A l3 8/8A a 8/8A v S. If we define 

p\{k - n - r 2 )p\{k + r 4 ) r fi (r 1 )S(k - r 1 )T v (r 2 )S(k - n - r 2 ) 

r a (r 3 )s(k + n)r^ n )S(k)\ (3.23) 

and use the result (|3.22| ) and the definitions (|3.13|) for the cutoff functions we find for the 
polarization tensor Fig.[3] 

6 See appendix C. 
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Figure 3: The representation of the amplitude A^ ua/3 (ri, r 2 , r 3 , r 4 ). 



r 5U( r i> r 2 , r 3 , r 4 ) = A fiua p(r 1 , r 2 , r 3 , r 4 ) + A^ain, r 2 , r 4 , r 3 ) + A^p^n, r 4 , r 2 , r 3 ) 

+A fiav p(r 1 ,r 3 , r 2 , r 4 ) + ^^(ri, r 3 , r 4 , r 2 ) 
+^ /1 ^(ri,r4,r 3 ,r 2 ). (3.24) 

Now we are ready to show that rjfj ~ is as expected finite 0] and is regular in the (IR) 



domain. After the substitution fl3.15| ) and a judicious shift of the integration variable in 



( |3.23j ), the product of the cutoff functions becomes proportional to p 6 {k), where p(k) is 
defined in (|3.12|) . From the rotational invariance of the integral and the properties (|A.1| ) and 
( [A.2j ) of the matrices, it follows that the asymptotic expression of ( |3.23| ) as — q 2 — > +oo 



is 



V«/3 ~ Tr (2r / ,r,r Q r^ - r^r^ - r r„r a v v ) J dk p^k)^^. (3.25) 

As a corollary it is readily seen that the potentially logarithmic divergent term of ( |3.24| ) 
vanishes identically in taking the sum, i.e. the polarization tensor is finite. Is the polarization 
tensor regular? 

The standard one loop expression of the polarization tensor aB which is known to 
be finite and gauge invariant |16| , can be deduced from fl3.23| ) and (|3.24j ) if one substitutes 



respectively the matrices by the Dirac's matrices 7^. As a result, the structure of the 



difference T^Lg — g is given by (|3.22|) after the substitution of the matrices T M by 



Dirac's matrices 7 M , if in this expression the last term in the right member is omitted. Since 
this difference is actually finite and gauge invariant, we can show along the same lines as 
before that it vanishes identically. Therefore the polarization tensor has the standard form, 
and thus is regular. Since for n > 2 the 2n-photon amplitudes are finite, we can show in a 
similar way that they are also regular. We end this section by discussing the status of the 
n-photon amplitudes when n is odd. 
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If S is the regularized electron propagator ( |3.6| ) and if the matrix T M is given by ( |2.23| ), 
the requirement of gauge invariance imposes that the n-photon amplitude, whose structure 
is similar to ( |3.22| ), must be proportional to the integral of the trace of n products of T^S. 
Under charge conjugation the matrices transform exactly in the same way as Dirac's 
matrices 7^. Hence it follows by virtue of Furry's theorem that the n-photon amplitude 
vanishes when n is odd. 



4 The mass operator and the vertex function 

First we derive the expression of the mass operator at the one loop level. 

For this sake we take the functional derivative of the equation of motion ( |2.33| ) defined 
by ( |2.20|) and ( |2.24j ) with respect to the C- number function ip b {p') and express the sources 
of the electron field according to the generating functional of (1PI) functions by relations 
( f2.31| ). We thus obtain the inverse of the full electron propagator ( |C.1| ) in terms of the 
mass operator £ as 

r£W) = (plW-m) ab S(p-p')-E ab (p,p') (4.1) 
Z ab (p,p') = -le |^Ti 2 & ) (A;^pOC(p,fc)^4^5 cd (A;,A ; ')e G 



+00 „ 

- E / dm ' dp " r 2V,p') F ^(p -^ k ~ k ')pl^ 

n=0 L 

+ (n + 2) J dq p 2 (k)p 2 (k - q) p 3 (q)T ^(q) SJ ^ q) F n+1 (p -k,k - k' - q) 

S cd (k', P ")e G + 0($. (4.2) 

Here the F n (p,q) which are defined in appendix B are expressed in terms of the functional 
derivative 5/5 J and the latin letters refer to the spin indices. The contribution of order 
e 2 which comes from the bracket of expression ( [4.2[) is proportional to T fl (q)q u D flu (q) and 
vanishes owing to the properties ( |A.2| ). Taking into account the fact that the functional 
derivative 5S/5J^ can be written in terms of the product of the full photon propagator 

D^(r,q) = l ^) \ J J & l { _ q) (4.3) 

times the functional derivative 5S/5A U whose expression at order e is given by ( |3.5|) , we get 
for the mass operator at order e 2 

E(p,p') = 5{p-p')i^ J dkD^(k-p)T,(p,k)S(k)T u (k,p). (4.4) 

This operator can be represented with the help of the vertex Fig.|I]a by the diagram of Fig.[|. 
In this expression D^ v (k) is the free photon propagator whose expression 

D^(k) = -^(v-pW-^)) ( 4 - 5 ) 
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P-k 



Figure 4: The mass operator. 



is deduced from fl2.21|) and T A1 (p, k) is defined in (|2.25|) . Owing to the property ( [2.26D of the 
matrix T^p, k), we can notice that we recover the standard non regularized form for the 
mass operator as the cutoff A tends to infinity. The mass operator E can be written as 



E(p)5(p-p') = E(p-p') 
S(p) = A(p 2 )+j)B{p 2 ). 



(4.6) 
(4.7) 



As stated before, we can replace freely the free propagator S{k) by the expression ( |3.15| ) 
without changing the result of integration in (|4.4j ). Then if we use definitions ( |3.13|) and 
( |3.12j ) for the cutoff functions and neglect finite terms of the order 0(oo/*) we get for the 
scalar functions A and B 



A{p 2 



p 2 B(p 2 



-im 



{2-kY 



dk \p 2 {k-p) 



[K 2 -(k-p) 2 Y 



+ 3p 2 (k-p)p(k) wz 



p 2 —2kp+m 2 



[K 2 ~{k-p) 2 Y{m 2 -k 2 ) 



(k-p) 2 ){m 2 ~k 2 ) 

(4.f 



l Wf J dkp 2 {k-p) {2pl(fe) [<ta _ (fc _ p ^ ](ma _ fca) + [kp(p2{k)-pi{k) 



+ 2p 2 



P 2(k) 



i 



[K 2 -{k-p) 2 Y 



(4.9) 



where the last integral of the right member of ((Of) is finite. In order to avoid the problem 
of potentially (IR) divergences due to the momentum carried by the internal photon line, 
we give an arbitrary small mass k to the photon. After Feynman's parametrization of the 
product of propagators, with the help of formulas ( p.6| ), ( |L).7| ) and ( p.9| ), the integration 
over the four momentum is easily performed and we obtain in the Feynman gauge the known 



results |l7j 

A(p 2 ) 



m- 



h, log ^ - 7 - 1 - \ log 54 - J dx log ■ 



-p x(l—x)+m (1—x)+xk 



B(P 2 ) = -f{log^-I-^-|log 



dx x log 



-p 2 x(l— x)+m 2 (l— x)+xk 2 
7fl 



(4.10) 



(4.11) 
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We now calculate the vertex function at the one loop level. If we define the (1PI) function 

r( 3 ) as, 



((27T 



,4\2 



5 3 r 



T^\p,p',r) = e8(p' - p-r)T^{p,p') 



(4.12) 



SA»(-r)84!(p)8ip(p') 

in the same way as we deduced the mass operator, we get from the equation of motion (|2.33| ), 
( CT and ( ggf 



rf (p,p',r) = 5(p' -p-ryT^p') + (2tt) z 



5Am (_ r ) 



;£(p,p') + £W))- (4-13) 



In this expression 

S(p,p') = ^ j dkdqD^{k-p A )Y a {p 1 k)S{k } p' + q)Yp{p' + q 1 p') (4.14) 



is the mass operator at the one loop order, T^p, p') is given by the definition (|2.25f) and 
E'{p,p') = J dkdk'dp" r$(pV)|e 2 \F 2 (p - k,k- k')p\{k)$ 



+ J dq p 2 (k)p 2 (k - q)p3(q)T v (q)F 1 (p - k, k - k' - q)-^ 



(-<?) 



S c d(k',p" 



+ e 



3 i 



\F z {p -k,k- k')p\{k)$ + Jdq p 2 {k)p 2 {k - q)p 3 (q)T u (q) 



F 2 (p — k,k — k' — q) 



8J v {-q) 



S cd (k',p") + 0(}*) f. 



(4.15) 



If we keep only the terms of order e 3 in the expression ( |4.15|) and use the recursion formulas 
( p.5|) we can show respectively that the first and the second bracket of £' vanish on the 
one hand after integration over dk and on the other hand after successive integration over 
dp",dk' and dk. Then according to the relation ( |3.5| ), the functional derivative of S ( [4.14| ) 
with respect to the C-number function N 1 is easily obtained. We thus get for the vertex 
function Fig.[5] at the one loop level 

ifW) = 7m + (§r J dk D^ik-p^^^S^T^k + ^Sik + r^fsik+r,/). (4.16) 



As the (UV) cutoff A tends to infinity due to the property ( |2.26|) , we recover formally 
the standard non regularized expression for the vertex function. 

Now we work in Feynman gauge and assume that the external electron lines are on 
mass-shell, thus we use freely the Gordon relation and express as usual the vertex function 
as 

rfW) = 7^(1 + ^1(0) + ^r u F 2 (r 2 ), (4.17) 
where r = p' — p is the momentum transfer. In this case the computation can be simplified 



if one notices that the term p\{p)^) — P\{k)^ which occurs in r^(p, k) ( |2.25|) is identical to 

S-\p) - S-\k) and that 



u(p)S-\p) = u(p)0{± 



qo \ 



(4.18) 
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Figure 5: The vertex function. 



p being the momentum of the outgoing electron. A similar reduction occurs for the term 
Fp(k+r, p')u(jp'). We give to the photon an arbitrary small mass k in order to avoid the (IR) 
divergence and substitute the free electron propagators by their expressions Q3.15|) . After 
Feynman parametrization of the product of propagators, we perform the integration over 
the internal momentum using the formulas of appendix C. Then, if we define the integrals 



h{r 2 ) 



I,(r 2 



tlx 



m 2 —r 2 x(l—x) 



l 

■ 2 x(\-x) 10 & 



h(r 2 ) 



m 2 —r 2 x(l—x) 



m 2 —r 2 x(l—x) 



(4.19) 



we get for the forms factors Fi(r 2 ) and F 2 (r 2 ) the expressions 



FAr 2 



a) 1 
7T I 2 



log £ - h - | - \ log 2 - \{2m 2 - r 2 )I 2 (r 2 ) + \ 



(2m 2 - r 2 ) log - 



+ (3m 2 



/i(r 2 ) - \h{r 2 ) 



(4.20) 

'2V ) = 2^" L J U' )■ ( 4 - 21 ) 
\i being an arbitrary unit of mass. As expected from ([4.17|) and ( [4.2 1| ) we recover the known 



F,(r 2 ) 



a -m 2 lAr 2 ) 



result for the anomalous magnetic moment of the electron |L7j . The expression of the form 
factor Fi(r 2 ) which is (IR) and (UV) divergent is also standard. 



What about the prediction of the (WT) identity? From the general relation ( 2.34 ) we 
can deduce the identity 

rT^ 3) (p,p', r) = ep 3 (r) r (2) (p + r, p ) — , p' — r) , r = p' — p. (4.22) 

If Z\ and Z 2 are respectively the renormalization constants of the vertex function and of the 
external electron lines associated with the mass operator i.e. 

7 _ 1 - §XM 



]/>=m 



(4.23) 
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the identity ( |4.22| ) implies that Z\Z^ X must be finite. If the external electron lines are 
off mass-shell, we can show that a term proportional to afr^/r 2 contributes to the vertex 
function too. Since the derivative of this term is not regular near the origin, the expression 
r^dY^' j dr p does not vanish as — ► 0. As a consequence, the (WT) (|4.22 ) does not imply 
the equality of Z\ and Z 2 . This explains the reason why Z\ differs from Z 2 in first order in 
a by a numerical constant as it can be easily seen. 



5 Conclusion 

We have constructed in four dimensions a continuous regularization scheme of QED based 
on a non local extension of the action. Since the measure of the path integral is invariant 
under the gauge group, this regularization scheme is actually non-perturbative. Once the 
regularized action is fixed by the choice of the cutoff functions, the regularized amplitudes 
can be calculated in a straightforward way from the path integral without the need to adjust 
some parameters or to define some integrals formally. In order to illustrate this fact we 
have deduced the regularized (1PI) functions at the one loop order from the equations of 
motion, which are now mathematically well-defined objects and thus can be represented by 
regularized Feyman diagrams. It follows from the expression Q2.12 ) of the C-number function 



Cn(z, z', x), which was imposed by the necessity of gauge invariance, that the theory keeps 
trace of the non-locality in some potentially divergent (1PI) functions through additional 
terms which are regular in the (IR) domain, but whose derivatives are not. This is the 
case for the polarization operator and for the vertex function if this latter is calculated off 
mass-shell. Although the polarization operator is transverse, it contains a quadratically 
divergent term whose derivative is not regular in the (IR) domain. Contrary to known gauge 
invariant regularization procedures, like the Pauli -Villars regularization where the quadratic 
divergent piece is not gauge invariant and then legislated to zero, this term must be removed 
by a specific choice of the counterterm S\. In the case of massless QED some terms of the 
non gauge invariant part of closed fermions loops are not unambiguously regularized to zero. 
This fact is reminiscent of the Nielsen-Ninomiya Theorem for chiral fermions on the lattice. 

Despite the fact that we recover the known result for the anomalous magnetic moment 
of the electron at the one loop order, we must investigate all the implications of the (IR) 
behavior of this regularized form of QED for physical process. These questions, the problem 
of non-perturbative renormalization, the possibility of dynamical mass generation for the 
photon through a new non-local interaction and the extension of this regularization scheme 
to non-abelian gauge theories like QCD, where non-perturbative effects are known to occur, 
will be investigated elsewhere. 
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Appendix A 

The matrices satisfy : 

{I»,r„(r)} = {I», Tl J = 2(^-^) (A.l) 

r M r M (r) = 0, l v T^r) lv = -2T,{r) (A.2) 

Appendix B 

The Fourier transform of e ieL ^ ,s/ ^ is 

F(p,q) = f dxdy e i{px+qy) e ieL{x ' y) . (B.l) 



If we define 

F n (p,q) = J dxdye l(px+ ^L n (x,y) (B.2) 



we have the expansion 



F M = E^nfo?) (B.3) 

n=0 

F (p,q) = ((2vr) 4 ) 2 5(p)5(q), (B.4) 
and the following recursion formulas hold for the F^s 

K+i(p, q) = -i J df p 3 (r)^ (F„(p -r,q)- F n (jp, q - r)) A»(r). (B.5) 
In addition one can easily show by induction that 

J dk k^F n (p- k,k-p') = 0, j dk F n (p - k,k-p') = 0, (B.6) 
for all p and p'. 

Appendix C 



If we take as usual the functional derivative of the electron source rj a (k) fl2.31p with respect 
to rjb(k'), the inverse of the full electron propagator (|3.3|) 



r^(k,k') = -(2*)' ,™ (c.i) 



5ipa{k)&i> b {k') 
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is defined by 

J dk" TW(k, k")S cb (k", k') = 6 ab 6(k - k'). (C.2) 
In the same manner the inverse of the full photon propagator (|4.3| ) 

r$(k,k') = (2n) \ A ^l { _ kl) (C.3) 

verifies 

J dk" r$(k, k")D au (k", k') = 6,J(k - k>). (C.4) 

The relations (|C.2|) and (|C.4|) are the basic tools for expressing, by means of Schwinger's 
technique of functional differentiation, general Green's functions in terms of (1PI) functions 
and propagators. For instance, taking the functional derivative of ( |C.2| ) with respect to the 
C- number function A^{— r'), we obtain the relation 

■ _ f dk'dk" SUp,k") 5r Z?"i? S da (k>,k). (C.5) 



Appendix D 

In this appendix we give the asymptotic form of the integrals which we encounter in the 
calculation of the relevant (1PI) functions. We start with the evaluation of the integral 

urh (-k 2 ~2qk+c) a ■ K^- 1 ) 

where e is an infinitesimal parameter. At first we make the shift of variable k — > k — q and 
we express the denominator of the integrand of ( p.l| ) by means of Schwinger's parametric 
integral. After a Wick rotation the integral over dk is easily performed ||, and we get 

T(a) J {ae+u) 2 V I 

Now if we expand exp [— e 2 (b — aq) 2 / (ae + «)] in power series and integrate over du we obtain 
the integral / in terms of the degenerate hypergeometric function ^(a, n, z) [Tj| as 



/ = m 2 e^ aq - 2bq \at) a - 2 |tf(a, a - 1, (q 2 + C)aej + ^-^t v {b ~ aq) 2lJ a~ v 

ty(a, a-l-v, (q 2 + C)ae^. (D.3) 
When n is an integer, the degenerate hypergeometric function ^(a, n, z) can be expanded 



as 



IS 



f n—l 

\B( a n z) = |y (-i) fc (fc-i)ir(a-fc) - k , y-na+ 

vyix, 10, 6) T(a-n+l)\ (n-k-l)\ F(a) Z ^ T(a 



r(a-n+l) I (n-fc-1)! r(a) 1 T(a) (n+k-l)\k\ 

k fc=l fc=0 

+ 1) +^(^ + fe) -tp(a + k) -logs]!, n>l (DA) 
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with the convention that for n = 1 the first sum is identically zero. In addition the relation 

m 

$>{a, n, z) = z l - n ^{a-n + l, 2 - n, z) (D.5) 



holds. In the above formulas, tp(x) is the ip function [18]. For a = 2, we use (|D.4|) and (|D . 5|) 
and the integral ( |D . 1| ) becomes 



dk 



a e(ak^+2bk) 



m z - log ae + - 1 - log(<f + C) + 0(e) log e 



(-k 2 -2qk+C) 2 

Similarly for the following integrals we obtain the final result 



dk A 



e e(ak A +2bk) 



(-k 2 -2qk+C) 2 



VK 



q^logae - \b^a 1 - qjip(l) - § - log(g 2 + C) 



-0(e) log e 



;d.6) 



(D.7) 



„ h h p ti(ak 2 +2bk) 
Ah. V" 6 

ar " (-k 2 -2qk+C) 2 



k k P c{ak 2 +2bk) 
Ah. "l'"" 6 

afh (-k 2 -2qk+C) 3 



VK 



iSV(ae) 1 - ^^(g 2 + C) + q^„J logae 
- aq)^b - aq) v oT 2 + \g^{q 2 + C) (V(l) - \ - log(g 2 + C) 
- aq) 2 a~ 2 + §((&- ag)^ + (6 - ag)^ ) a -1 



|^( logae - V>(1) + § + log(g 2 + C) 



(D.8) 



f4^ + 0(e)loge 



(D.9) 
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